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Abstract 

We study the quantum correlation and classical correlation dynamics in a spin-boson model. For 
two different forms of spectral density, we obtain analytical results and show that the evolutions 

Qh| of both correlations depend closely on the form of the initial state. At the end of evolution, all 

-I— ' ■ 

C ' correlations initially stored in the spin system transfer to reservoirs. It is found that for a large 






family of initial states, quantum correlation remains equal to the classical correlation during the 
course of evolution. In addition, there is no increase in the correlations during the course of 



1> ' evolution. 
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It is generally accepted that entanglement plays a crucial role in quantum information 
processing. However, it is fragile due to the inevitable interaction of the system with its 



a. 



environment. There have been many works on entanglement for all kinds of systems 
Many extraordinary results have been obtained, such as entanglement sudden death (ESD) 
4i^, and sudden birth [7|]. 

However, more and more works show that, even without entanglement, there are quantum 
tasks |8l-ll0| superior to their classical counterparts. This is due to their nonzero quantum 
correlation (Q). Recently, some special results have been obtained for Q and its counterpart 
classical correlation (C), such as sudden change ll|, |l2|. Therefore, it is important and 
interesting, to investigate Q and C in more open systems. Q, introduced by Henderson and 
Vedral [l3|, is identical to quantum discord introduced by Olliver and Zurek \i^ in the 



case of two qubits ll|, ll5|, ll6|. Q describes all nonclassical correlations in a state, while 
entanglement is only a special part of the correlations. 



To obtain Q, one must acquire C first. As shown in |17|, C{psj^s2) = max{nj}[5'(psj — 
>S'{n }(Psi|s2)]; where the maximum is taken over the set of positive-operator-valued mea- 
surements (POVM) {Hj} in partition S2. Q is then obtained by subtracting C from the 
quantum mutual information. In fact, Q is asymmetric 1^, it depends on which partition 
the measurement is carried out in. Only under the condition that S{pa) = S{pb) (where 
A and B are the two partitions, and S{p) denotes the Von Neumann entropy of state p), 
does it become symmetric [l7J. It is difficult to obtain Q for a general state, owing to the 
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18| on the dynamics of correlations. 



maximal operation. There are only a few works |ll . 

and more work is needed to gain a deep understanding of quantum and classical correlations. 

In this Brief Report, we investigate the dynamics of Q and C in a spin-boson model, 
which is conventionally used to study open system. For two different forms of spectral 
density, we acquire analytical results for both Q and C for different partitions. It is shown 
that evolutions of Q and C depend closely on the form of the initial state. For a large family 
of initial states, Q remains equal to C during the course of evolution. We also find oscillation 
of Q and C. Furthermore, we find that all Q and C initially stored in the spin system run 
into the boson reservoirs continuously without sudden death, which is in sharp contrast 
with the case for entanglement. In addition, the absence of an increase in correlations is 
investigated. 

The model we are considering is two spins interacting independently with their own boson 



reservoirs without interaction between the two spins [19]. The Hamiltonian reads as 

2 N 



i=l k=l 

2 N 



(1) 



j=l fc=l 



(with A^ — )■ oo). Here, bik {i = 1,2) is an annihilation operator of the fcth model in 
the zth reservoir with corresponding frequency Uk- cr^, a^ (o'j^), and u represent the Pauh 
operator, raising (lowering) operator, and Zeeman splitting of the ith spin, and Qi^k denotes 
the coupling strength between the fcth model in the ith reservoir and the corresponding 
spin. As is well known, all information on the coupling between the system and reservoir is 
included in the spectral density, written as J{u)) = J2k \9k\'^^i.^ " ^k) (Where we suppose 
both spins and reservoirs are identical). 

Let us first consider the evolution of the initial state with two excitations in the spin 
system, 

|<l>o) = (a|0),jO),, + /3|l),Jl),J|0),jO),„ (2) 

with |ap + |/3p = 1, where the collective state \Q)ri = Y[k=i lOfc)^ indicates there is no 
excitation in reservoir Tj. Evolution of the whole system is given by 

N 



k=l 

N 



X (e(t)|i)s.|o)., + 5^AKt)|o),,|i,).,; 



1=1 



+ a|0),jO),,|0),jO),,. 



(3) 



After introducing the collective state 11)^. = - X]a;=i '^fc(^)|lfe)r-i with one excitation in the 
reservoir, the reduced density matrix of the spin system is 

/|ap+|/3xT a/3*e^\ 

1/3^x1' 

1/3^x1' 

V a*pe |/3eT/ 

In our case, we find that the maximum can be obtained at ^ = 7r/4 ^]; therefore, the 



PsisM 



(4) 



analytical results of C and Q are 



C{Ps.sAt)) = g(p,,,,(t)) = H(|/3eP) 



1 



- H(-(l + v/l-4|/3exP)), (5) 

where we have introduced the Shannon entropy H(a;) = — a;log2a; — (1 — a;)log2(l — x). For 
the reservoirs, C and Q read as 

CiPr,r,{t)) = Q(p.,.,(t))=H(|/3x|') 



- H(-(l + v/l-4|/3exP)). (6) 

First, when the spectral density is fiat, 

Jiu) = 7, (7) 

where 7 is a constant that is commonly used as the Markov approximation with the interval 
of the spectral density much broader than the corresponding energy scale of the system. We 
obtain parameters 

X = VT^^^, e = e-^*/^ (8) 

The concurrence l20| of the system is max{0, e^^*(|a;/3| — |/3p(l — e^"^*))}. If |a| < |/3|, 
there will be ESD [19']. As time approaches infinity, the asymptotic behavior of C for the 
spin system p^^^j becomes 

hm C(p,,,,(t)) ~ (1/3^ - |/3r)7te-^^7ln2, (9) 

which does not depend on the initial value C(psiS2(0)) = —(1 — l/5niog2(l ~ l/^P) ^ 
|/3plog2|/3p. It is shown that C of the spin system is depleted gradually without sudden 
death. 

According to Eq. (5), C is equal to Q during the evolution. As shown in Figure 1(a), 
with a = l/v2, /3 = l/v2, and a = 1/vlO, /3 = 3/vlO, the correlations tend to zero 
monotonically and continuously. 

According to Eq. (6) for reservoirs, C equals Q, and their asymptotical behaviors are 
described by 

lim Qiprrr,it)) = Cipr,r,it)) ~ Qips.sM)- (lO) 

Figures 1(a) and (b) show that both Q and C initially stored in the spin system run 
into reservoirs gradually as the whole system achieves equilibrium. Evolutions of Q and 
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FIG. 1: (Color online) Evolution of Q and C among different partitions with flat spectral density for 
the initial state (a|0)sj0)s2 + /?|l)si|l)s2)|0)ri|0)r2- Blue diamonds and magenta squares denote 
Q and C for the Bell state' respectively, while dark triangles and red circles denote Q and C 
for a = l/vTo, /3 = S/VTo. (a) Spin system ps^s^^- (b) Reservoirs Pr^^r-i- (c) Spin s\ with its 
corresponding reservoir ri, Psin- (d) Spin s\ with the other reservoir r2, Ps^T-i,- 




FIG. 2: (Color online) Dynamics of Q and C for the initial state (a|0)s^ |0)s2 +/3|l)si |l)s2)|0)ri |0)r2 
and the spectral density taking the Lorentz form W jX = \/200 . Blue diamonds and magenta 
squares denote Q and C for the Bell state, respectively, while dark triangles and red circles denote 
Q and C for a = l/\/IO, /? = 3/-v/l0- (a) Spin system PsiS2- (b) Reservoirs Prir2- (c) Spin si with 
its corresponding reservoir ri, Ps^n- (d) Spin si with the other reservoir r2, Psir2- 



C between the system and reservoirs may reflect the transference processes of Q and C 
between the system and reservoirs. Si (^2) becomes correlated with ri (r2) owing to their 
interaction as shown in Fig. 1(c), therefore, ri tangles with r2 through Si and 82- At the end 
of evolution there is no entanglement between reservoirs and system; therefore, the system 



can be discarded without any effect, and all correlations transfer to reservoirs. 

We also find that Q, C, concurrence between si amd ri, and concurrence between S2 
and r2, do not increase during the course of evolution. They are measured by the square 
of correlations Q{tf = Q{Ps,s,{t)f + Q{ps,r,{t)f + Q{ps,rMf + Qipr.r^it)? < Q(0)2, 

C{tf = Cips.sM? + CiPs.rM? + Cip,,rAt)? + C{pr,rM? < C^(0)^ and 0071(1)' = 

Con{ps^s2if)Y + Con{ps^r2if)Y + Con{ps^ri{t)Y + C on{pr^r2if)Y ^ Con(O)^. It can be seen 
that for partitions si and S2, interaction between S2 and r2 is simply a local unitary operation; 
therefore, correlations will not exceed those initially stored between si and S2- 

Until now, we have mainly focused on the weak coupling regime, where dynamics of 
correlations are monotonic and irreversible. Things will be different in the strong coupling 
regime. We consider spectral density taking the Lorentz form, which is usually encountered 



in a cavity 2l|. Here we take the spectral density centered at loq (resonance with the spin 
system) : 

where 1/A is used to measure the correlation time of the reservoir. For A ~ 0, we obtain 
J{uj) ~ W''5{uj — Wo), which displays a strong non-Markov effect. We have the parameters 



^ = e-^*/2[A/VA2 - 4Py2sinh(VA2 - AWH/2) 



+ cosh(VA2-4Vr2^/2)], 

X = yr^. (12) 

Concurrence due to the strong non-Markov effect (measured by W/X) is given by 
max{0, |q;*/3^^| — |/3^^^X^|}) which has oscillation. C and Q oscillate too, as shown in Fig. 
2(a). According to Eq. (5), for this kind of initial state, although the coupling has changed, 
C is equal to Q. As shown in Fig. 2(a), the evolutions of Q and C are continuous. Evo- 
lutions of Q and C for reservoirs also have oscillation, as shown in Fig. 2(b). From the 
analytical results, we find that both Q and C initially stored in the spin system transfer to 
the reservoirs over time. Figures 2(c) and 2(d) show the evolutions of correlations for Psin 
and Psir2: respectively. A non increase in correlations and concurrence is noted, but both 
correlations and concurrence oscillate. 

We then consider the initial state with only one excitation in the spin system: 

|$o) = (a|0),Jl),,+/3|l),jO),J|0),jO),,. (13) 
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FIG. 3: (Color online) Evolution of Q and C for the initial state (a|0)sjl)s2 +/3|l)si|0)s2)|0)ri|0)r2 
with both reservoirs having a flat spectrum. Blue diamonds and magenta squares denote Q and 
C for the Bell state, respectively, while dark triangles and red circles denote Q and C for a = 
l/-v/TO, /3 = 3/vTO. (a) Spin system PsiS2- (b) Reservoirs />rir2- (c) Spin si with reservoir ri, Ps^r^- 
(d) Spin si with reservoir r2, Psxr2- 



(a) 


-•-Q, r<-1/2"' 




HD-C, a = 1/2"' 


' 11 


^ii^Q,a = l/10"^ 


lill 


-^C,a=1/10"' 



2.0 2.5 3.0 







(d) ^•^Q,aM/2"^ 




[1 * i^ ^C,a=1/2"= 




1 


-A-Q, a=1/10"^ 








-•-C, a=1/10'^ 


CD 

i-0.2 

o 
o 

0.1 


Hi 


1 


Ima 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



FIG. 4: (Color online) Dynamics of Q and C for the initial state (a|0)sjl)s2 +/3|l)si|0)s2)|0)rj0)r2 
with the spectral density of reservoirs taking the Lorentz form W/\ = -v/200. Blue diamonds and 
magenta squares denote Q and C for the Bell state, respectively, while dark triangles and red circles 
denote Q and C for a = l/\/lO, /3 = 3/\/l0- (a) Evolution of spin system Ps-^S2- (t>) Reservoirs 
Pr-ir2- (c) Spin si with reservoir ri, ps^rx- (d) Spin si with reservoir r2, Psir2- 



Evolution of the spin system is given by the reduced density matrix: 

l\x? o\ 



Ps^S2{t) 



(14) 



la^P «/3iel' 

a*m^ m? 

\ 0/ 

Interestingly, in this situation, C has the same value as that in the first situation given 
by Eq. (5), while for Q, we have 



QiPs.s^t)) = -H(|eP) + H(|<r) 

+ H(i(l + v/l-4|/3exP)). 

For reservoirs Prir2) we obtain 

Cipr.rM) = H(|/3xr)-H(i(l-Vl-4|/3exP)), 

QiPr^rM) = H(^(l + v/l-4|/3exP))-H(|xP) 
+ H(|axl'). 



(15) 



(16) 



With fiat spectral density, we obtain the parameters in Eq. (8). Concurrence of the spin 
system reads as \a* (3\e~'^'^ , and therefore, there is no ESD. As time approaches infinity, Q 
can be represented asymptotically by 



limg(p,,,,(t))~g(0)e 



-7* 



(17) 



withg(O) =H(|a|2). 

For this kind of initial state, Q is no longer equal to C as in Fig. 3(a), but the correlations 
fall to zero monotonically and continuously. For reservoirs, both correlations approach a 
definite value monotonically as in Fig. 3(b). The transference processes of quantum and 
classical correlations are similar to those of the initial state with two excitations, as shown 
in Figs. 3(c) and 3(d). As time tends to infinity, both Q and C initially stored in the spin 
system run into reservoirs. In addition, there are no increases in the correlations. 

For the Lorentz form in Eq. (11), we have the same parameters as in Eq. (12). Concur- 
rence of the spin system is given by |q;*/3^^|. Q is not equal to C. As shown in Fig. 4(a), 
correlations of the spin system also oscillate. Figure 4(b) shows evolutions of correlations in 



reservoirs. Figures 4(c) and 4(d) show evolutions for siVi and Sir2. Although there is large 
oscillation due to the non-Markov effect, there are no increases in correlations. 

In conclusion, we have investigated the dynamics of both Q and C in the spin-boson model 
with two spins independently coupled to their own reservoirs. From the analytical results 
obtained for Q and C, we studied the evolutions of Q and C among different partitions in 
detail. We found that the dynamics of Q and C depend closely on the form of the initial 
state. At the end of evolution, all Q and C initially stored in the spin system transfer to 
reservoirs. We found that Q is more robust than entanglement, for there is no sudden death 
with Q. At the same time, during the evolution process, all partitions had nonzero Q, which 
is not the case for entanglement. We also found that for a large family of initial states, Q 
remains equal to C during the course of evolution. In addition, there was no increase in 



either correlation. There are many differences between this paper and Ul, ll2|]; for example, 
we focused on a more concrete case of the spin-boson model, which is usually encountered, 
and we concentrated on the transference processes of both Q and C. 
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